Resumen

Row

Última actualización de datos realizada el 25 de febrero de 2021 con base en los datos publicados por la Dirección de de Agua del Ministerio de Ambiente y Energía de Costa Rica.

Row

código del punto muestreado

PS-ESQ-01

Cuerpo de Agua

Piedras Blancas

Cuenca

Esquinas

Coordenadas

583319,21

Coordenadas

970837,56

Row

Provincia

Puntarenas

Cantón

Osa

Distrito

Piedras Blancas

Row

Resumen de resultados de los índices de la calidad del agua superficial

Utilizando el Índice Holandes y el Índice Biológico (BMWP-CR), según el Decreto Ejecutivo N°33903-MINAE-S
Mes Año Régimen_de_lluvia Indice_Holandés Interpretacion_IH Indice_Biológico Interpretacion_IB
Agosto 2019 Estacion lluviosa 4 Contaminacion incipiente 187 Aguas de calidad excelente
Noviembre 2019 Transición Lluviosa-Seca N/A N/A N/A N/A
Febrero 2020 Estación Seca N/A N/A N/A N/A
Mayo 2020 Transición Seca-Lluviosa N/A N/A N/A N/A

Row

Mapa de ubicación del punto de muestreo

Información de los datos publicados

1. El monitoreo de los cuerpos de agua superficiales y la clasificación de su calidad se realizó según:

  • Decreto N° 33903-MINAE-S. REGLAMENTO PARA LA EVALUACIÓN Y CLASIFICACION DE LA CALIDAD DE CUERPOS DE AGUA SUPERFICIALES.

  • Plan nacional de monitoreo de la calidad de los cuerpos de agua superficiales.

  • Programa nacional de monitoreo de la calidad de los cuerpos de agua del país.

2. Para una mejor interpretación de los resultados, se aclara:

Según el artículo 20 del Decreto N° 33903-MINAE-S, en aquellos casos en que el valor del Índice de Clasificación Holandés y el valor del Índice Biológico presenten una diferencia de más de una clase entre ellos, el muestreo se deberá repetir por una vez, en un período no mayor a 30 días naturales. En caso de que se mantenga la diferencia de clases se clasifica en la de menor calidad.

Row

Notas:

  1. En los casos donde el resultado se reporta comio N/A, es debido a que no se tienen resultados del Índice por la imposibilidad de tomar la muestra, ya sea por el acceso al sitio de muestreo o porque no discurría agua por el cauce del cuerpo de agua.

  2. Para mayor información sobre el monitoreo de la calidad de los cuerpos de agua superficiales, escríbanos al correo:

CONSULTA DE DATOS

Row

CONSULTA DE DATOS DEL MONITOREO EN COSTA RICA (2015-2020)

Column

Datos del monitoreo entre los años 2015 al 2020

Distribución espacial de registros

---
title: "MONITOREO DE LA CALIDAD DEL AGUA SUPERFICIAL (ETAPA I) DE COSTA RICA"
output: 
  flexdashboard::flex_dashboard:
    orientation: rows
    social: menu
    source_code: embed
    vertical_layout: scroll    
---

```{r setup, include=FALSE}

#-------------------- Paquetes --------------------

library(flexdashboard)
library(plotly)
library(dplyr)
library(tidyr)
library(sf)
library(leaflet)
library(rgdal)
library(knitr)



#-------------------- Colores ---------------------

color_camp_1 <- 'blue1'
color_camp_2 <- 'blue'
color_camp_3 <- 'blue'
color_camp_4 <- 'blue'

color_info <- '#9ccf75'
color_ubic <- '#84c057'


#--------------------- Íconos ---------------------

icono_camp_1 <- 'fas fa-hand-holding-water'
icono_camp_2 <- 'fas fa-fill'
icono_camp_3 <- 'fas fa-flask'
icono_camp_4 <- 'fas fa-clipboard-list'
icono_info1 <- 'fas fa-cogs'
icono_info2 <- 'fas fa-cogs'
icono_info3 <- 'fas fa-edit'
icono_info4 <- 'fas fa-clipboard-check'

#--------------- Otros parámetros -----------------

# Separador para lectura de datos CSV
caracter_separador <- ','
```


```{r, include=FALSE}
#--------------- Archivos de datos ----------------

archivo_datos_PNM <- read.csv('https://raw.githubusercontent.com/EsmeVar/PNM_FASE_5/main/2021_04_28_PNM_INFO_F5.csv', sep = ";")

#convertí aparte la tabla en un shape, es una opción en QGIS y la agrego directamenete y le doy la proyección
PNM <- st_read("C:/Users/usuario/Documents/ESMERALDA/Año 2021/PNMI/PUBLICACION_DATOS_PNM_FASE_5/DATOS_PNM/PNM_WGS84_2.shp") %>%
  st_transform(4326) %>%
  st_make_valid()

#para crear un mapa sería bonito agregar las cuencas y hacer coropletas
cuencas <- st_read("C:/Users/usuario/Documents/ESMERALDA/Año 2021/PNMI/PUBLICACION_DATOS_PNM_FASE_5/DATOS_PNM/Cuencas hidrograficas/CUencas hidrograficas.shp")%>%
  st_transform(4326) %>%
  st_make_valid()

names(archivo_datos_PNM)[1] <- "nombre" #cambia el nombre de la columna cuencas para que coincida con la del shp de cuencas

#buscar y reemplazar para que coincidan los nombres entre ambas bases 
archivo_datos_PNM$nombre <- gsub("Baru", "Barú", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Chirripo", "Chirripó", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Curen/a", "Cureña", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Frio", "Río Frío", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Grande de Terraba", "Térraba", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Jesus Maria", "Jesús María", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Moin", "Moín", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Peninsula de Nicoya", "Península de Nicoya", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Peninsula de Osa", "Península de Osa", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Reventazon", "Reventazón", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Sarapiqui", "Sarapiquí", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Tarcoles", "Tárcoles", archivo_datos_PNM$nombre)
archivo_datos_PNM$nombre <- gsub("Estrella", "La Estrella", archivo_datos_PNM$nombre)

#interseca los puntos con el shape para mezclar la información
sf_datos <- inner_join(cuencas, archivo_datos_PNM,  by = "nombre")

```


Resumen
=======================================================================
Row {data-height=10}
-----------------------------------------------------------------------
### **Última actualización de datos realizada el 25 de febrero de 2021 con base en los [datos publicados por la Dirección de de Agua del Ministerio de Ambiente y Energía de Costa Rica](http://mapas.da.go.cr/mapnew.php)**.




Column {.sidebar}
-----------------------------------------------------------------------


[![Dirección de Agua](C:/Users/usuario/Documents/ESMERALDA/Año 2021/PNMI/PUBLICACION_DATOS_PNM_FASE_5/LOGOS DA-MINAE/da-logo230x.png)](http://www.da.go.cr) [![MINAE](C:/Users/usuario/Documents/ESMERALDA/Año 2021/PNMI/PUBLICACION_DATOS_PNM_FASE_5/LOGOS DA-MINAE/Logo MINAE 2020_230x.jpg)](https://minae.go.cr/)























**Elaboración:** Esmeralda Vargas Madrigal **Colaboración:** Kevin Row {data-height=130} ----------------------------------------------------------------------- ### código del punto muestreado {.value-box} ```{r} valueBox(value = paste(format(archivo_datos_PNM[1,3])), caption = "Código del sitio de muestreo", icon = icono_info1, color = color_info ) ``` ### Cuerpo de Agua {.value-box} ```{r} valueBox(value = paste(format(archivo_datos_PNM[1,2])), caption = "Cuerpo de Agua", color = color_info ) ``` ### Cuenca {.value-box} ```{r} valueBox(value = paste(format(archivo_datos_PNM[1,1])), caption = "Cuenca", color = color_info ) ``` ### Coordenadas {.value-box} ```{r} valueBox(value = paste(format(archivo_datos_PNM[1,12])), caption = "Longitud (Coordenadas CRTM 05)", icon = icono_info3, color = color_info ) ``` ### Coordenadas {.value-box} ```{r} valueBox(value = paste(format(archivo_datos_PNM[1,13])), caption = "Latitud (Coordenadas CRTM 05)", icon = icono_info3, color = color_info ) ``` Row ----------------------------------------------------------------------- ### Provincia {.value-box} ```{r} valueBox(value = paste(format(archivo_datos_PNM[1,15])), caption = "Provincia", color = color_ubic ) ``` ### Cantón {.value-box} ```{r} valueBox(value = paste(format(archivo_datos_PNM[1,16])), caption = "Cantón", color = color_ubic ) ``` ### Distrito {.value-box} ```{r} valueBox(value = paste(format(archivo_datos_PNM[1,17])), caption = "Distrito", color = color_ubic ) ``` Row {data-height=210} ----------------------------------------------------------------------- ### **Resumen de resultados de los índices de la calidad del agua superficial** ```{r} knitr::kable( archivo_datos_PNM[1:4, 18:24 ], caption = '**Utilizando el Índice Holandes y el Índice Biológico (BMWP-CR), según el Decreto Ejecutivo N°33903-MINAE-S**' ) ``` Row {data-width=400} ----------------------------------------------------------------------- ### **Mapa de ubicación del punto de muestreo** ```{r} mapa1 <- leaflet() %>% addTiles() %>% addProviderTiles(providers$Stamen.TonerLite, group = "Stamen Toner Lite") %>% addProviderTiles(providers$Esri.WorldImagery, group = "Imágenes de ESRI") %>% addMarkers(data= PNM, group = "PNM", popup= paste( "Sitio de muestreo: ", PNM$Codigo )) %>% addLayersControl( baseGroups = c("OpenStreetMap", "Stamen Toner Lite", "Imágenes de ESRI") ) # Despliegue del mapa mapa1 ``` ### **Información de los datos publicados** **1. El monitoreo de los cuerpos de agua superficiales y la clasificación de su calidad se realizó según:** - Decreto N° 33903-MINAE-S. REGLAMENTO PARA LA EVALUACIÓN Y CLASIFICACION DE LA CALIDAD DE CUERPOS DE AGUA SUPERFICIALES. - Plan nacional de monitoreo de la calidad de los cuerpos de agua superficiales. - Programa nacional de monitoreo de la calidad de los cuerpos de agua del país.

**2. Para una mejor interpretación de los resultados, se aclara:** Según el artículo 20 del Decreto N° 33903-MINAE-S, en aquellos casos en que el valor del Índice de Clasificación Holandés y el valor del Índice Biológico presenten una diferencia de más de una clase entre ellos, el muestreo se deberá repetir por una vez, en un período no mayor a 30 días naturales. En caso de que se mantenga la diferencia de clases se clasifica en la de menor calidad.

Row {data-width=50} ----------------------------------------------------------------------- ### **Notas:** 1. En los casos donde el resultado se reporta comio N/A, es debido a que no se tienen resultados del Índice por la imposibilidad de tomar la muestra, ya sea por el acceso al sitio de muestreo o porque no discurría agua por el cauce del cuerpo de agua. 2. Para mayor información sobre el monitoreo de la calidad de los cuerpos de agua superficiales, escríbanos al correo: daguas@da.go.cr CONSULTA DE DATOS ======================================================================= Row {data-height=20} ----------------------------------------------------------------------- ### CONSULTA DE DATOS DEL MONITOREO EN COSTA RICA (2015-2020) Column {data-height=1000} ----------------------------------------------------------------------- ### Datos del monitoreo entre los años 2015 al 2020 ```{r} sf_datos %>% st_drop_geometry() %>% select(Cuenca = nombre, Código = Codigo, "Cuerpo de agua" = Cuerpo.de.Agua, Provincia = Provincia, Cantón = Canton, Distrito = Distrito) %>% DT::datatable(rownames = FALSE, options = list(searchHighlight = TRUE, language = list(url = '//cdn.datatables.net/plug-ins/1.10.11/i18n/Spanish.json') )) ``` ### Distribución espacial de registros ```{r} leaflet() %>% fitBounds(lng1 = -86, lng2 = -82, lat1 = 8, lat2 = 11) %>% addProviderTiles(providers$OpenStreetMap.Mapnik, group = "Open StreeT Map") %>% addTiles(urlTemplate ="https://mts1.google.com/vt/lyrs=s&hl=en&src=app&x={x}&y={y}&z={z}&s=G", attribution = 'Google', group = "Google Maps") %>% addMarkers(data= PNM, group = "PNM", popup = paste( "Sitio de muestreo: ", PNM$Sitio.de.M, "
", "Código: ", PNM$Codigo, "
", "Cuerpo de agua: ", PNM$Cuerpo.de, "
", "Cuenca: ", PNM$Cuenca, "
", "Fase de implementación: ", PNM$Fase, "
", "Año de muestreo inicial: ", PNM$An.O.INICI, "
", "Año de muestreo final: ", PNM$An.O.FINAL, "
", "Provincia: ", PNM$Provincia, "
", "Cantón: ", PNM$Canton, "
", "Distrito: ", PNM$Distrito, "
")) %>% addPolygons( data = sf_datos, stroke=T, fillOpacity = 0, color="red", weight=0.8, opacity= 2.0, group = "Cuencas", popup = paste( "Cuenca: ", sf_datos$nombre, "
", "Provincia: ", sf_datos$Provincia, "
", "Cantón: ", sf_datos$Canton, "
", "Distrito: ", sf_datos$Distrito, "
")) %>% addLayersControl(baseGroups = c("Google Maps", "OpenStreetMap"), overlayGroups = c("PNM", "Cuencas"), options = layersControlOptions(collapsed = TRUE)) %>% addScaleBar() %>% hideGroup(c("Cuencas")) %>% addMiniMap( toggleDisplay = TRUE, position = "bottomleft", tiles = providers$OpenStreetMap.Mapnik) ```